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I. 


Introduction 


When  a  metal  liner  is  driven  to  a  velocity  V  by  an  explosive  charge, 
it  is  often  important  to  know  not  only  the  magnitude  of  its  velocity 
but  also  the  direction  in  which  each  liner  element  moves.  This  direc¬ 
tion  is  conveniently  defined  by  the  angle  6  between  the  velocity  vector 
of  the  liner  element  and  the  perpendicular  to  the  initial  liner  surface. 
We  shall  call  6  the  liner  projection  angle.  For  a  fragmentation  charge 
(or  exploding  bomb),  the  projection  angle  will  determine  the  final  angu¬ 
lar  distribution  of  the  scattered  fragments.  For  explosive  cladding  or 
shaped-charge  liner  collapse,  it  will  determine  the  collapse  angle  6 
of  the  given  charge.  The  specified  qualities  of  either  the  cladding 
bond  or  the  jet,  crucially  depend  on  this  collapse  angle. 

A  formula  for  determining  the  projection  angle  6  from  the  velocity 
of  the  liner  V  and  the  velocity  U  by  which  the  detonation  wave  front 
sweeps  past  the  liner  surface  was  first  proposed  by  Taylor  [1]: 

sin5  =  V/2U  (1) 

This  formula,  which  has  been  the  most  extensively  used  one  to 
date  (see  References  2-8),  is,  however,  accurate  only  under  steady  state 
conditions  where  the  detonation  wave  sweeps  past  identical  cross  sec¬ 
tions  of  the  explosive-liner  geometry.  For  non-steady  cases,  the  Taylor 
formula  is  not  applicable  since  the  existence  of  either  a  gradient, 

9V  Bt 

—  ,  of  the  velocity  V  along  the  liner  or  a  gradient,  ,  of  the  typical 

acceleration  duration  t(Z  denotes  the  length  tangential  to  the  liner 
surface)  may  very  significantly  affect  the  angle  6,  sometimes  causing 
very  big  deviations  between  the  Taylor  predictions  and  experimental 
measurements . 

A  more  accurate  formula  for  the  liner  projection  angle  is  there¬ 
fore  needed.  It  will  lead  to  a  more  accurate  description  and  under¬ 
standing  of  processes  such  as  the  collapse  of  the  conventional  shaped 
charge  or  the  hemispherical  liners  and  the  formation  of  the  self-forg¬ 
ing  fragment. 

Recently,  Randers-Pehrson  [9]  derived  a  non-steady  liner  projection 
angle  formula  by  empirically  f-ftting  a  formula  to  the  numerically  cal¬ 
culated  results. 

In  the  present  paper  an  analytically  derived  formula  is  obtained. 

It  is  based  on  the  assumptions  that  (1)  the  detonation  pressure  acts 
normally  on  the  liner,  and  (2)  the  angle  (0-6),  which  will  be  defined 
later,  is  small. 

We  found  that  for  the  small  angle  assumption  to  be  valid,  the 
initial  radius  of  curvature  must  be  small  in  comparison  with  the  dis¬ 
tance  traveled  by  the  liner  during  acceleration.  A  third  assumption  is 
that  the  internal  forces  in  the  liner  metal  can  be  neglected.  The  new 
formula  was  tested  against  both  numerical  calculations  and  experimental 
data  and  was  found  to  predict  the  projection  angle  6  more  accurately 
than  the  steady  Taylor  formula. 
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In  this  report  we  address  only  the  equation  governing  the  angle  6 
and  the  velocity  parameters  V0  and  t.  No  attempt  is  made  here  to 
determine  values  of  VQ  and  t  as  functions  of  explosive  and  liner  geometry 
and  properties.  It  is  obvious  that  values  of  Vq  and  x  are  needed  for 
the  eventual  application  of  this  equation.  A  simple  method  in  deter¬ 
mining  Vq  and  x  will  be  the  objective  of  future  research.  For  the 
present  purpose  of  ascertaining  the  accuracy  of  the  unsteady  Taylor 
angle  equation,  we  shall  use  values  of  Vq  and  x  determined  either  by 
two-dimensional  computer  code,  or  by  experimental  measurement. 


I I .  Derivation  of  Basic  Equations 

We  shall  first  derive  the  differential  equations  describing  the 
liner's  motion.  We  shall  limit  ourselves  to  liners  with  very  large 
initial  radii  of  curvature.  In  other  words,  the  liner  has  an  almost 
straight  formation  line.  The  liner  could  either  be  an  axisymmetric 
shell,  or  a  "plane  strain"  plate.  The  explosion  wave  front  is  assumed 
to  be  cylindrically  symmetric  for  the  shell  case  with  its  axis  of  sym¬ 
metry  coinciding  with  that  of  the  liner's.  We  denote  by  &  the  length 
coordinate  along  the  liner  and  by  U  the  velocity  by  which  the  detona¬ 
tion  wave  front  sweeps  past  the  liner  surface  along  the  direction  of 
the  liner  surface. 

The  first  differential  equation  deals  with  the  increase  of  liner 
velocity  with  time.  Let  V  and  6  be  the  magnitude  and  direction  of  the 
liner  velocity  at  the  point  Z  on  the  liner  at  a  specific  time,  t.  We 
make  the  assumption  that  the  gas  pressure  always  acts  perpendicular 
to  the  liner  surface.  Whan  a  liner  segment  elongates  or  shrinks  as  it 
is  pushed  or  pulled  by  its  neighboring  segments,  a  force  component  in 
the  liner  direction  may  also  appear.  We  assume  that  this  force  can  be 
neglected  during  the  acceleration  time.  We  denote  by  9  the  angle  be¬ 
tween  the  original  liner  formation  line  direction  and  the  current  for¬ 
mation  line  direction,  at  SL.  Then,  we  can  see  from  Fig.  1  that  the 
additional  velocity  vector  d$,  at  a  given  position  l ,  induced  by  the 
force  acting  perpendicularly  to  the  liner  during  the  time  dt,  adds  to 
the  current  velocity  vector  ^(t)  to  form  the  new  velocity  vector 
^(t+dt) .  As  a  result,  the  velocity  increases  in  magnitude  by  the 
amount  dV  =  d | V |  =  |dv|  cos (0-6)  where  |d^|  is  the  magnitude  of  the 
above  mentioned  additional  velocity  vector.  At  the  same  time,  the 
tangential  velocity  component  |d^|  sin(9-6)  causes  the  velocity  direc¬ 
tion  to  change  by  the  an(,ie  d6 .  We  can  relate  the  angle  d6  to  this 
component  by  the  equation: 

|  dV J  sin(0-6)  =»  V*d6. 

Substituting  |dv|  =  dV/cos(Q-6)  we  get  the  scalar  equation  for 
the  velocity  magnitude: 

dV  tan(9-6)  «*  V*d 6  .  (2) 

Dividing  by  dt  and  denoting  the  differentiation  with  respect  to  the 
time  by  a  dot  above  the  symbol  we  finally  obtain 

•  v 

6  »  —  tan (9-6) .  (3) 

The  next  equation  describes  the  influence  of  the  existence  of 
a  velocity  gradient  along  the  liner.  We  make  the  assumption  that  a 
liner  will  not  elongate  while  being  accelerated,  or  that  its  elonga¬ 
tion  may  be  neglected. 

We  are  interested  in  calculating  the  rate-of-change  of  liner 
39 

slope  .  Referring  to  Fig.  2,  we  note  that  during  the  time  dt, 
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Initial  Liner  Contour 


Liner  a 
vclocit 
^nd  cha 
V. 


point  A  will  move  the  distance  V(t)dt  to  point  C  in  a  direction  making 
the  angle  9(H)  -  5(H)  with  the  normal  to  the  current  liner  surface 
at  point  A.  Point  B  will  move  at  the  same  time  the  distance 

dH^  •  dt  to  point  D.  The  angle  at  A  (see  Fig.  2)  between  the 
initial  and  current  liner  direction  is  equal  by  definition  to  9(H). 

The  angle  d9  between  CD  and  AB  is  determined  by  the  components 
of  V  dt  and  (V  +  dV)dt  normal  to  the  current  liner  at  H  thus: 


d9  m  tan(d9) 


(4) 


At  the  limit  of  small  dH  this  yields: 

9  ■  -  V'  cos(9-6)  (5) 

where  the  prime  denotes  differentiation  with  respect  to  H,  and  dot 
with  respect  to  time. 


HI.  Integration  of  Basic  Equations 

Eqs.  (3)  and  (5)  represent  two  equations  governing  the  three 
quantities  V,  9,  and  6.  Assuming  that  the  spatial  and  temporal  dis¬ 
tribution  of  V  for  a  given  problem  is  known,  we  can  solve  for  0  and  6. 
Under  the  approximation  tan(9-6)  *  9-6.  Eq.  (3)  now  becomes 

•  • 

•  v  V 

s  +  |-«-7e-  <6> 

Eq.  (6)  is  a  first  order  differential  equation  which  for  the 
initial  conditions  0=6=V=O  at  t  <_  T  has  the  general  solution: 


6  = 


ft 


9V  dt 


0-^1  V0  dt 


(7) 


where  T  is  the  instant  when  the  liner  begins  to  accelerate.  From 
Eq.  (5)  we  obtain  under  the  above  approximation: 


0  =  - 


V' cos (0-6) dt  *  -  V'dt. 


T 


(8a) 


When  we  substitute  0  from  Eq.  (5),  the  second  term  of  Eq.  (7)  becomes 
V0dt  =  ^  |  V  V'  cos(0-6)dt  =  ^  |  W  dt.  (8b) 


We,  therefore,  obtain  the  general  solution  for  6 
6  -  -  |  V'dt  +  ~  |  (V2)'dt. 


(9) 


N 


IV.  Application  to  Exponential  Acceleration 

In  order  to  apply  Eq.  (9)  to  a  practical  problem,  we  shall  assume 
for  V  the  exponential  form: 


V(£,t) 


v0(t) 


1  -  exp 


[-(*]} 


when  Vq,T,  and  r  are  functions  of  £  only.  Though  this 
fit  perfectly  for  all  cases,  it  is  reasonably  accurate 
have  studied  and  easily  integrable  for  the  calculation 
Eq.  (9).  The  corresponding  acceleration  is 


(10) 

form  may  not 
for  cases  we 
of  6  given  by 


3V  v0  r  /  t-T\~] 

a,7F“TexpL"  Wj 

The  quantity  T(£)  is  the  arrival  time  of  the  detonation  wave  at  the 
point  £  on  the  liner,  t  is  the  characteristic  acceleration  time  of 
this  point  and  Vq(£)  is  the  final  asymptotic  velocity  reached  by  the 
liner.  The  time  t-0  is  taken  when  the  detonation  wave  reaches  a  con¬ 
venient  point  on  the  metal  say  £-0. 


To  facilitate  the  integration,  we  interchange  the  order  of  time 
integration  and  differentiation  with  respect  to  £. 


6 


V  dt  + 


11^ 

2V  3£  J 


2 

V  dt. 


(ID 


(the  differentiation  with  respect  to  T  vanishes  since  V  ■  0  at  t  ■  T) . 
For  the  V  function  given  by  Eq.  (10), 

/■t 

V  dt  -  VQ(t-T)  -  tV  (12) 


/: 


v2dt  -  VQ[V0(t-T)  -  TV]  -  J  TV2. 


(13) 


Substituting  Eq.  (12)  and  (13)  into  Eq.  (11)  yields: 


S  -  £{TV-V0(t-T)]  -  £  {V0CTV-V0(t-T»  +  i  TV2)  . 


(14) 


To  compare  this  equation  with  experimental  results,  we  shall  be 
Interested  in  the  value  6  will  obtain  as  t  -*■  “.  At  this  limit;  V  be¬ 
comes  equal  to  Vq,  V'  becomes  equal  to  Vq'  and  the  expression  for  6 
becomes: 
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m 


mMaummaam 


as 


V  T' 

6  -  - 4  TV’  +  7 

2  2  0  4  0 


(15) 


Let  U  be  the  detonation  wave  sweep  speed,  then  T'  is  equal  to 
1/U  for  a  detonation  on  the  aris  of  symmetry  and  Eq.  (15)  becomes: 


-0  -  I  TV  ’ 

2U  2  0 


+  i 

+  4  1 


u’ 


(16) 


The  first  term  is,  for  a  small  angle  6,  identical  to  the  Taylor 
formula,  Eq.  (1),  and  the  other  two  are  the  correction  terms. 
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V.  Constant  Acceleration  Equation 


In  deriving  Eq.  (16)  we  assumed  the  exponential  fora  for  the 
velocity  time  dependence.  Even  though  this  form  very  well  simulates 
the  experimental  measurements  it  is  not  the  only  form  one  can  choose 
to  fit  the  experimental  data.  When  solving  Eq.  (9)  with  a  constant 
acceleration  assumption  at  the  period  of  time  T  <  t  _<(T  +  T^)i.e., 


*  - 


we  find 


and  for  the  form: 


V.T '  .  . 

_ A  _  v  »  +  A  T  *  v 

2  6  cv0  6  c  0 


,  1/2 


we  obtain 


V  ,  , 

_0 _ A  T  V  *  +  —  t'  V 

2  6  SR  V0  12  SR  0 


(17) 


(18) 


(19) 


(20) 


It  can  easily  be  seen  that  the  igR  parameter  of  Eq.  (19)  ex¬ 
actly  corresponds  to  3t,  when  t  is  the  corresponding  parameter  in 
the  exponential  acceleration  (Eq.  10). 

Using  an  acceleration  fora  such  as  Eq.  (17)  or  Eq.  (19)  may  be 
more  convenient  in  obtaining  a  closed  form  solution  in  more  compli¬ 
cated  cases,  e.g.,  when  the  effect  of  the  liner  curvature  on  the 
angle  6  cannot  be  neglected. 
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VI .  Comparison  with  Previous  Work  and  Two-Dimensional  Code  Calculations 


Let  us  now  compare  Eq.  (16)  with  the  equation  obtained  empiri¬ 
cally  by  Randers-Pehrson  [9]  using  Two-D  code  calculations  and 
experiment.  The  formula  given  in  Ref.  9  is  as  follows: 


6 


—  -  —  tV  '  -  —  (tV  ’)^. 
2U  2  0  5  ^  0  ’ 


(21) 


(Note  that  Eq.  (21)  agrees  with  Eq.  (16)  in  the  first  two  terms). 

To  perform  this  comparison  we  need  to  know  Vq  and  Vq',t  and  t' . 
There  exists  some  experimental  data  for  Vq  and  Vq'  from  exploding 
cylinder  tests.  Experimental  data  for  t  (except  in  special  steady 
state  cases)  seems  lacking,  however.  We,  therefore,  performed  Two-D 
code  simulations  to  calculate  Vq,Vq',t,t'  and  6  for  two  specific 
geometries.  Our  general  approach  is  to  use  output  from  a  Two-D  code 
of  Vq  and  t  at  a  number  of  £  locations  to  numerically  calculate  Vq' 
and  t'.  Then  substituting  these  values  of  Vq  and  Vq',  t  and  x ' 
into  Eq.  (16),  we  obtain  a  value  for  6  which  we  shall  denote 
(F  -  formula).  The  value  of  <5p  is  then  compared  with  the  Two-D  code 
value  of  6  denoted  Sc(c-code).  Wherever  experimental  measurements  of 
6  are  available  for  a  specific  charge  we  shall  denote  them  6ex.  We 
expect  to  find  a  good  agreement  between  6C  and  6ex  when  using  a 
reliable  Two-D  code.  Therefore,  a  check  against  experimental  data  is, 
in  fact,  a  check  of  the  Two-D  code  validating  its  use  for  checking 
Eq.  (16).  For  completeness,  these  values  are  also  substituted  into 
Eq.  (21)  and  the  result  denoted  6Rp 

The  two  dimensional  code  we  found  most  convenient  to  apply  to 
this  purpose  is  a  Lagrangian  code  named  TEMPS  in  which  the  explosive 
is  treated  by  a  two-dimensional  finite-difference  grid  similar  to 
HEMP,  TOODY  and  other  Two-D  Lagrangian  codes.  The  liner  in  TEMPS, 
however,  is  described  as  an  array  of  mass  points  having  tensile  and 
bending  forces  between  each  two  neighboring  points.  The  liner  is 
thus  being  treated  as  one-dimensional.  This  not  only  saves  compu¬ 
tation  time  but  also  facilitates  the  data  reduction  from  the  code  out¬ 
put. 

A.  Conical  Shaped  Charge 

The  first  configuration  studied  is  the  81.3  mm  diameter,  42°, 
conical  shaped  charge  depicted  in  Fig.  3.  In  this  example  we  choose 
£=0  at  the  cone  liner  apex  and  £  is  the  distance  from  this  point  along 
the  formation  line.  The  charge  is  initiated  by  plane  detonation  at 
the  rear  of  the  explosive. 

From  the  TEMPS  code  calculations,  the  velocity  history  for  each 
liner  mass-point  position  £  is  recorded  and  plotted.  A  typical  plot 
is  shown  in  Fig.  4.  The  asymptotic  final  velocity  reached  is  defined 
as  Vq(£)  for  each  mass  point.  These  are  tabulated  in  Table  1  for 
this  case.  This  final  velocity  is  then  plotted  as  a  function  of  £ 
as  shown  in  Fig.  5a.  This  curve  is  then  numerically  differentiated 
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Aluminum  Casing 


Lgure  3.  Drawing  of  81.3mm  diameter  42°  conical  shaped  charge. 


using  the  simple  difference  formula  AVq/A2  to  obtain  Vq' (Z)  which  is 
also  contained  in  Table  1. 

Next  we  must  estimate  the  value  of  t  for  each  mass  point.  This 
was  obtained  by  fitting  Eq.  (10)  to  the  TEMPS  velocity  data.  The 
arrival  time  T  used  in  Eq.  (10)  is  the  theoretical  arrival  time  at  the 
explosive  metal  interface  and  obtained  by  dividing  the  distance  be¬ 
tween  the  plane  of  initiation  and  the  point  Z  by  the  detonation 
velocity  Up.  As  is  indicated  in  Fig.  A,  the  wave  arrival  time  pre¬ 
dicted  by  the  TEMPS  code  doe9  not  coincide  with  the  calculated  theo¬ 
retical  arrival  time.  This  inaccuracy  in  the  code  simulrtion  is 
caused  by  the  code's  smearing  of  the  detonation  front  in  the  finite- 
difference  calculation  scheme.  Therefore,  in  performing  the  fit  for  t, 
we  did  not  choose  x  as  the  value  of  t  when  V  =  Vf)(l-e~l)  =  0.632  Vq, 
as  is  indicated  by  Eq.  (10).  This  value  of  x  would  be  inaccurate 
since  it  is  close  to  and  heavily  influenced  by  the  code's  calculation 
of  wave  arrival  time.  We  found  it  more  reliable  to  fit  the  curves 
closely  in  the  region  V  .86  Vq.  The  values  of  t  obtained  from  this 
procedure  are  given  in  Table  1.  These  values  of  x  are  also  plotted 
as  a  function  of  Z  in  Fig.  5b.  A  line  is  fitted  through  the  points 
and  x'  calculated  by  simple  numerical  differentiation.  The  values 
are  shown  in  Table  1.  We  see  that  for  this  charge  x  is  approximately 
constant  and  therefore  x'  is  a  very  small  quantity. 

These  values  are  then  substituted  into  Eq.  (16)  and  Eq.  (21)  as 
well  as  Taylor's  relation  and  compared  to  the  code  calculations.  The 
quantity  6C  is  the  angle  of  the  velocity  vector  as  given  by  TEMPS 
when  the  liner  reaches  a  final  velocity  Vq.  The  detailed  data  are 
given  in  Table  1  and  plotted  graphically  in  Fig.  6.  We  see  that  the 
Eq.  (16)  gives  a  significant  improvement  over  the  Taylor  relation 
and  agrees  well  with  the  TEMPS  code.  Note  that  since  the  last  term  in 
both  Eqs .  (16)  and  (21)  is  small,  the  difference  between  Eq.  (16/  and 
Eq.  (21)  is  small  for  this  example. 

B.  Exploding  Cylinder 

The  second  configuration  studied  is  an  exploding  cylinder  which 
is  identical  lo  the  example  used  by  Randers-Pehrson  [9,10].  It  con¬ 
sists  of  a  steel  pipe  segment  101.6  mm  in  length  and  a  diameter  of 
50.8  mm  filled  with  0CT0L  as  shown  in  Fig.  7.  For  this  example, 
direct  experimental  measurements  of  6  and  Vq  are  presented  in  Ref.  10. 
In  these  experiments,  the  expanding  cylindrical  pipe  breaks  into  frag¬ 
ments.  The  speed  and  direction  of  motion  of  the  fragments  are  measured 
by  means  of  x-ray  shadowgraphs  taken  at  several  predefined  times. 
Similar  experiments  were  also  conducted  in  References  11  and  12. 

To  provide  a  complete  verification  of  Eq.  (16),  however,  the 
quantity  x  is  also  needed.  Since  no  acceleration  data  was  measured 
in  these  experiments,  we  cannot  provide  an  experimental  value  for  x 
and  hence  a  full  experimental  verification  of  Eq.  (16)  at  this  point 
in  time  Is  not  possible.  In  fact,  as  mentioned  in  the  previous  sec¬ 
tion,  measurements  of  metal  motion  during  acceleration  have  been  re¬ 
ported  for  only  a  very  special  steady  state  case  (see  [13,14)  in 
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POSITION  ALONG  LINER  l (mm) 

Figure  5b.  Plot  showing  t  as  a  function  of  location  l  from  TEMPS  calculat 
conical  charge. 
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which  Vg*  and  t'  are  identically  zero.  Therefore*  in  our  present 
comparison,  TEMPS  code  calculations  will  again  be  used  to  supply  the 
quantities  t  and  x'  (as  well  as  Vg  and  Vg*).  The  experimental  data 
will  thus  be  used  to  verify  the  final  values  for  Vg  and  6  given  by 
the  code  calculations  as  well  as  the  6  given  by  Eqs.  (IS)  and  (21). 

The  results  of  the  TEMPS  calculations  for  Vg  and  6  are  compared 
to  the  experimental  data  in  Figs.  8  and  9.  Also  shown  in  these  fig¬ 
ures  are  the  HEMP  simulation  results  presented  by  Randers-Pehrson  in 
[9].  Both  TEMPS  and  HEMP  results  are  within  the  spread  of  the  experi¬ 
mental  data. 

The  procedure  to  determine  x,  x',  Vq  and  Vq*  from  the  TEMPS  data 
was  described  in  the  previous  section.  A  velocity  vs.  time  plot  for 
a  typical  mass  point  is  shown  in  Fig.  10.  A  plot  of  x  vs.  £  is  given 
in  Fig.  11.  The  values  of  all  the  necessary  quantities  are  listed  in 
Table  2.  These  values  were  then  substituted  into  Eqs.  ^16Tiand  (21) 
to  compute  and  Sop.  Fig.  12  is  a  comparison  of  <5p,  6gp,  the 
Taylor  formula  and  tne  TEMPS  result  6C.  We  see  that  considerable 
improvement  in  the  estimation  of  the  angle  6  is  achieved.  Again*  we 
note  that  the  third  term  in  both  Eqs.  (16)  and  (21)  is  small,  there¬ 
fore,  6p  and  6jjp  do  not  differ  significantly. 


POSITION  a (mm) 

Comparison  of  TEMPS,  HEMP  and  experimental  data  for  V_  of  an  exploding  cylinder 


exp Iodine  cylinder. 


O  temps 
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Figure  10.  Exponential  fit  of  TEMPS  Mass  Point  velocity  V  vs.  tine  at  £,=  -64 . 38mm. 


lot  showing  t  as  a  function  of  location  l  fitted  from  TEMPS  calculation  of 
:xploding  cylinder  charge. 
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Figure  12.  Comparison  of  projection  angled  calculated  from  formulas  and  TEMPS  Code 


TABLE  2.  RESULTS  OF  EXPLODING  CYLINDER  CALCULATIONS 
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VII.  Summary  and  Conclusions 

In  this  report  a  formula  for  the  explosive-metal  Taylor  angle  6 
in  the  unsteady  case  was  derived  from  basic  physical  principles  under 
the  assumptions 

(a)  the  explosive  r  ressure  always  acts  normal 
to  the  current  liner  surface. 

(b)  total  angular  motion  of  the  metal  is  small 
(i.e.,  9-6  is  small) 

(c)  forces  in  the  liner  are  small  during  the 
acceleration  time. 

Also,  the  present  case  is  restricted  to  liners  whose  meridional 
curvature  is  small.  Under  the  assumption  of  an  exponential  decaying 
acceleration  in  time  of  the  metal,  this  formula  is  given  by 


6 


V 


0 


where  Vg  is  the  final  velocity,  t  the  characteristic  acceleration  time, 
and  U  the  detonation  sweep  speed;  the  prime  indicates  spatial  dif¬ 
ferentiation  along  the  liner.  We  point  out  that  the  first  term  repre¬ 
sents  the  steady-state  Taylor  result.  The  remaining  two  terms  repre¬ 
sent  the  unsteady  effects.  Note  that  the  first  two  terms  are 
identical  to  the  semi-empirical  formula  given  in  [9]. 

for  the  collapse  of  a  conical  shaped  charge  and  the  explosion  of 
a  metal  cylinder  we  have  found  this  formula  to  yield  accurate  results 
when  compared  to  Two-D  hydrocode  calculations  and  to  provide  a  signifi¬ 
cant  improvement  over  the  Taylor  relation.  The  experimental  data 
available  also  shows  that  the  formula  is  accurate. 

More  recently,  liners  with  large  amounts  of  curvature  in  the 
meridional  direction  (formation  line)  such  as  hemispheres  have  shown 
promise  as  candidates  for  certain  warhead  applications.  For  such 
cases,  the  present  formula  is  not  applicable.  It  is  therefore  suggested 
that  the  current  analysis  be  extended  to  include  these  more  complex 
geometries . 
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1  Commander 

Naval  Research  Laboratory 
Washington,  DC  20375 

1  USAF/AFRDDA 

Washington,  DC  20311 

1  AFSC/SDW 

Andrews  AFB 
Washington,  DC  20311 

1  US  Air  Force  Academy 

ATTN:  Code  FJS-41  (NC) 

Tech  Lib 

Colorado  Springs,  CO  80840 

1  AFATL/DUR  (J.  Foster) 

Eglin  AFB,  FL  32542 

1  AFWL  (SUL,  LT  Tennant) 

Kirtland  AFB,  NM  87116 


1  AFLCy^WWMC 

Wright-Patterson  AFB,  OH  45433 
1  Director 

Lawrence  Livermore  Laboratory 
ATTN:  Dr.  J.  Kury 
P.  0.  Box  808 
Livermore,  CA  94550 

1  Director 

Lawrence  Livermore  Laboratory 
ATTN:  Dr.  M.  Wilkins 
P.  0.  Box  808 
Livermore,  CA  94550 

1  Director 

Lawrence  Livermore  Laboratory 
ATTN:  Dr.  E.  Lee 
P.  0.  Box  808 
Livermore,  CA  94550 
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No.  of 

Copies  Organization 

1  Director 

Lawrence  Livermore  Laboratory 
ATTN:  Dr.  H.  Horning 
P.  0,  Box  808 
Livermore,  CA  94550 

1  Director 

Lawrence  Livermore  Laboratory 
ATTN:  Dr.  J.  Knowles 
P.  0.  Box  808 
Livermore,  CA  94550 

1  Director 

Lawrence  Livermore  Laboratory 
ATTN:  Dr.  M.  Van  Thiel 
P.  0.  Box  808 
Livermore,  CA  94550 

1  Director 

Lawrence  Livermore  Laboratory 
ATTN:  Dr.  C.  Godfrey 
P.  0.  Box  808 
Livermore,  CA  94550 

1  Director 

Lawrence  Livermore  Laboratory 
ATTN:  Tech  Lib 
P.  0.  Box  808 
Livermore,  CA  94550 

1  Battel le-Columbus  Laboratories 
ATTN:  Mr.  Joseph  E.  Backofen 
505  King  Avenue 

Columbus,  OH  43201 

2  Dyna  East  Corporation 
ATTN:  P.  C.  Chou 

J.  Car leone 
227  Hemlock  Road 
Wynnewood,  PA  19096 


No.  of 

Copies  Organization 

1  Firestone  Defense  Research 

and  Products 

Division  of  the  Firestone 
Tire  and  Rubber  Company 
ATTN:  R.  Berus 
1200  Firestone  Parkway 
Akron,  OH  44317 

2  Honeywell,  Inc. 

Government  and  Aeronautical 

Products  Division 
ATTN:  C.  R.  Hargreaves 
R.  S.  Kensinger 
600  Second  Street 
Hopkins,  MN  55343 

2  Sandia  Laboratories 

ATTN:  Dr.  W.  Herrman 
Dr.  J.  Asay 

Albuquerque,  NM  87115 

1  Shock  Hydrodynamics 

ATTN:  Dr.  L.  Zernow 
4710-4716  Vineland  Avenue 
North  Hollywood,  CA  91602 

1  Systems,  Science  5  Software 

ATTN:  Dr.  R.  Sedgwick 
P.  0.  Box  1620 
La  Jolla,  CA  92037 

4  University  of  California 

Los  Alamos  Scientific  Lab 
ATTN:  Dr.  J.  Walsh 
Dr.  R.  Karpp 
Dr  C.  Mautz 
Tech  Lib 
P.  0.  Box  1663 
Los  Alamos,  NM  87545 

1  University  of  Denver 

Denver  Research  Institute 
ATTN:  Mr.  R.  F.  Recht 
2390  S.  Univeasity  Blvd. 
Denver,  CO  80210 
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No.  of 

Copies  Organization 

2  University  of  Illinois 
Dept  of  Aeronautical  and 
Astronautical  Engineering 
ATTN:  Prof.  A.  R.  Zak 
Prof.  S.  M.  Yen 
Urbana,  IL  61801 

Aberdeen  Proving  Ground 

Dir,  USAMSAA 
ATTN:  DRXSY-D 
B.  Oeheil 
G.  Johnson 
DRXSY-MP,  H.  Cohen 
Cdr,  USATECOM 

ATTN:  DRSTE-TO-F 
Dir,  USACSL 

Bldg.  E3S16,  EA 

ATTN:  DRDAR-CLB-PA 
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USER  EVALUATION  OF  REPORT 


{’lease  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number^ _ _ _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.)  _  _  _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.)_ _ _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name : 


Telephone  Number: 
Organization  Address: 


